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Scaling models of random N x N hermitian matrices and passing to the limit N oo leads to 
integral operators whose Fredholm determinants describe the statistics of the spacing of the eigenval- 
ues of hermitian matrices of large order. For the Gaussian Unitary Ensemble, and for many others 
as well, the kernel one obtains by scaling in the "bulk" of the spectrum is the "sine kernel" . 
"q^" ■ Rescaling the GUE at the "edge" of the spectrum leads to the kernel Ai'(ij^-Ai'(3:) Ai(ij) ^^^^.^ 

^\ , is the Airy function. In previous work we found several analogies between properties of this "Airy 

0^ ' kernel" and known properties of the sine kernel: a system of partial differential equations associated 

with the logarithmic differential of the Fredholm determinant when the underlying domain is a union 
(-H , of intervals; a representation of the Fredholm determinant in terms of a Painleve transcendent in 

' the case of a single interval; and, also in this case, asymptotic expansions for these determinants 

and related quantities, achieved with the help of a differential operator which commutes with the 

■ integral operator. In this paper we show that there are completely analogous properties for a class of 
' kernels which arise when one rescales the Laguerre or Jacobi ensembles at the edge of the spectrum, 

namely 

(N 

■ where Jaiz) is the Bessel function of order a. In the cases a = these become, after a variable 
nJ ' change, the kernels which arise when taking scaling limits in the bulk of the spectrum for the 

^ , Gaussian orthogonal and symplectic ensembles. In particular, an asymptotic expansion we derive 

■ will generalize ones found by Dyson for the Fredholm determinants of these kernels. 



Oh! 



I. INTRODUCTION AND STATEMENT OF RESULTS 
A. Introduction 



Scaling models of random N x N hermitian matrices and passing to the limit — > oo leads to integral operators 
• whose Fredholm determinants describe the statistics of the spacing of the eigenvalues of hermitian matrices of large 
^ ' order [ p^|j2^ . Which integral operators (or, more precisely, which kernels of integral operators) result depends on the 
matrix model one starts with and at which location in the spectrum the scaling takes place. 



For the simplest model, the Gaussian Unitary Ensemble (GUE), and for many others as well (see, e.g., 1 16 17| j2i^ , p4| ) , 
the kernel one obtains by scaling in the "bulk" of the spectrum is the "sine kernel" 

sin7r(x — y) 
tt{x - y) 

Precisely, this comes about as follows. If {0fe(x)}^g is the sequence obtained by orthonormalizing the sequence 
{x^ e~'^^/'^} over (— c!0,oo) and if 

N-l 

KN{x,y) = ^ (l)k{x) (j)k{y) , (1-1) 

fe=0 



*e-mail address: tracy@itd.ucdavis.edu 
^e-mail address; widom@cats.ucsc.edu 



1 



then in the GUE the probabihty density that n of the eigenvalues (irrespective of order) he in infinitesimal intervals 
about xi, . . . ,Xn is equal to 

Rn{xi,...,Xn) = det {KN{Xi,Xj)) ij=i^,,,^ri- 

The density of eigenvalues at a fixed point z is Ri{z), and this is ~ \/2N /n as N oo. Rescaling at z leads to the 
sine kernel because of the relation 

n TTX , Try sin7r(x-y) 
lim Ky(7 H , ^ + ^ — 



Rescaling the GUE at the "edge" of the spectrum, however, leads to a different kernel. The edge corresponds to 
z ~ V2N, at which point the density is ^ 2^Nb ^ and we have there the scaling limit P,0,E2| 



1 , / X I y , Ai(x)Ai'(w) - Ai'fx) Ai(w) 

N^oo22Ns 22 Nb 22 Nb x-y 

where Ai is the Airy function. In previous work |28j we found several analogies between properties of this "Airy kernel" 
and known properties of the sine kernel: a system of partial differential equations associated with the logarithmic 
differential of the Fredholm determinant when the underlying domain is a union of intervals p5[: a representation of 
the Fredholm determinant in terms of a Painleve transcendent in the case of a single interval |15|; and, also in this 
case, asymptotic expansions for these determinants and related quantities achieved with the help of a 

differential operator which commutes with the integral operator. (See |p7|| for further discussion of these properties 
of the sine kernel.) 

In this paper we show that there are completely analogous properties for a class of kernels which arise when one 
rescales the Laguerre or Jacobi ense mble s at the edge of the spectrum. For the Laguerre ensemble the analogue of 
the sequence of functions {0fc(a;)} in (1.1) is obtained by orthonormalizing the sequence 

over (0, cx)) (here a > —1), whereas for Jacobi one orthonormalizes 

{x''{l-xT'^{l + xf'''} 

over (—1,1). (Here a,/3 > —1.) In the Laguerre ensemble of (positive) hermitian N x N matrices the eigenvalue 
density satisfies [0,p3i, for a fixed x < 1, 



This limiting law is to be contrasted with the well-known Wigner semi-circle law in the GUE. The new feature here 
is the "hard edge" for x ~ 0. At this edge we have the scaling limit ||ll|] : 

where Ja{z) is the Bessel function of order a. Both limits follow from the asymptotic formulas for the generalized 
Laguerre polynomials. (Scaling in the bulk will just lead to the sine kernel and scaling at the "soft edge," x ~ 1, will 
lead to the Airy kernel.) The same kernel arises when scaling the Jacobi ensemble at —1 or 1. (Recall that in the 
Jacobi ensemble both ±1 are hard edges; see e.g. 23|.) 

For later convenience we introduce now a parameter A and define our "Bessel kernel" by 

K{x,y):=\ ^- ^— ^- — , x^ty) (1.2a) 

2[x-y) 

= ^(j„(\/i)2- J„+i(VJ)J„_i(V^)) (x = y). (1.2b) 
Before stating our results, we mention that in the cases a — we have, when A = 1, 

2^yKix\y^) = '^^^±^^, (1.3) 
7r(x — y) 7T(x + y) 

which are kernels which arise when taking scaling limits in the bulk of the spectrum for the Gaussian orthogonal and 
symplectic ensembles [Q. In particular, an asymptotic expansion we derive will generalize ones found by Dyson |^ 
for the Fredholm determinants of these kernels. 
We now state the results we have obtained. 
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B. The System of Partial Differential Equations 

We set 

J:= |J(a2j_i,a2j) K > 0) (1.4) 

and write D{J; A) for the Fredholm determinant of K (the operator with kernel K{x, y)) acting on J. If we think of 
this as a function of a = (ai, . . . , a2m) then 

2m 

dlog D{J; A) = - ^hiy R{aj,aj) daj (1.5) 
i=i 

where R{x, y) is the kernel of K{I — K)^^. We introduce the notations 

4>{x) VA J„(\/J), V(a;) a;0'(2;) (1.6) 

and the quantities 

: = (/ - A^-V (a,), := (/ - i^)" V (a,), (j = 1, . • . , 2m) (1.7) 

= (0,(/-i^)-V), (<^,(/-if)-V), (1.8) 
where the inner products refer to the domain J. The differential equations are 

(.^fc) (1.9) 

oak a-j - ak 

^^^_,^kM!li^p, (,^fc) (1.10) 

oak a-j - ak 



daj 4 ^ flj-afe 



du 



daj 

dv 
daj 



i-iyq^ (1.13) 
i-iypjqj. (1.14) 



Moreover the quantities R{aj,aj) appearing in (1.5) are given by 

a, R{a„a,) = ^(-1)^= a. ^"^P' -J^'^'^' + - \ [a' - a, + 2v)q^ + \ p, u. (1.15) 
fc#j aj ak 4 z 

These equations are quite similar to eqs. (1.4)-(1.9) of pS], as is their derivation. 



C. The ordinary differential equation 

For the special case J = (0, s) the above equations can be used to show that q{s]\), the quantity q of the last 
section corresponding to the endpoint s, satisfies 

s{q^-\){sq'y ^ q{sq'Y + \{s^a^)q+-^sq^{q^-2) ('= ^). (1.16) 

with boundary condition 
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g(s; A) 



2"r(l + a) 



(1.17) 



This equation is reducible to a special case of the Py differential equation;^ explicitly, if q{s) = (1 + y{x))/{l — y{x)) 
with s — x"^, then y{x) satisfies Py with a' = — = 7' = and S' — —2. (We have primed the usual Py 

parameters to avoid confusion with the a in our kernel. We mention that this special Py can be expressed algebraically 
in ter ms o f a third Painleve transcendent and its first derivative jl^ . We mention also that an argument can be given 
that (1.16) must be reducible to one of the 50 canonical types of differential equations found by Painleve, without an 
explicit ver ificat ion being necessary. This will be discussed at the end of section II B.) It is sometimes convenient to 
transform (I.IC) by making the substitution 

q{s) — cosTp{s), 



so that tp satisfies 



s 



1 

8^ 



sin (2^) 



cos -0 
4s2 sin^-f/) 



The Fredholm determinant is expressible in terms of q by the formula 

D{J; A) = exp 



log - qitf dt 



(1.18) 



(1.19) 



Denoting by R{s) minus the logarithmic derivative of D{J; A) with respect to s, we have also the representation 

2 _2 



1 



R[s) = - cos"' V(s) + « ^ IT ^(^ 



ds 



4s 



(1.20) 



Furthermore, R{s) itself satisfies a differential equation which in the Jimbo-Miwa-Okamoto a notation for Painleve 
III (see, in particular, (3.13) in |1J]) is 



[sG"f + (t' ((T - sa') (4a' - I) - [a'f = 



(1.21) 



where <j{s) = sR{s)] it has small s expansion 



where 



ct(s;A) = CaS^+" 
1 



1 - 



3 + 2a 



2(2 + a) 16(l + a)(2 + a)(3 + a) 



,2 2+2a 



1 + a 
1 

"(TT"a)2-" 



3 + 2a 
1 : tttS 



2(2 + a) 



(1.22) 



22"+2r(l + a)r(2 



We mention that in the special case a = and A = 1 we have D{J; 1) = e q{s; 1) = 1, ip{s, 1) = 0, and 
a(s,l) = s/4 exactly @,|,|lll. 



^The Painleve V differential equation is 

dx^ \2y y — 1 J \dxj x dx x'^ \ y J x 3/^1 

where a', /?', 7', and S' are constants. 
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D. Asymptotics 



Again we take J = (0, s) and consider asymptotics as s — > oo. From the random matrix point of view the interesting 
quantities are 



E{n;s) 



1)" (9" 



D{J-\) 



(1.23) 



This is the probabiUty that exactly n eigenvalues lie in J . The asymptotics of £^(0, s) = D{J; 1) are obtained from 
(1.19) using the asymptotics of q{s;l) or equivalently a{s;l) obtained from (1.16) or (1.21), respectively. (Our 



derivation is heuristic since as far as we are aware the corresponding Painleve connection problem has not been 
rigorously solved.) We find that as s — > oo, 



E{0;s) 



-s/4+a^ 



a _ 



(1.24) 



9a^ n , 3a 

s" + ( 

128 428 



51a^ 
1024 



1275a4 



1024 32768 



where Tq is a constant which cannot be determined from the asymptotics of q (or a) alone. However, as we mentioned 
above, when a = this expansion must agree with those obtained from formulas (12.2.6) of ||l^ (see also (12.6.17)- 
(12.6.19) in ijl^) after replacing s by n^t^. This leads to the conjecture 

_ G(l + a) 
~ (2^)"/2 

where G is the Barnes G-functio n pj . This conjecture is further supported by numerical work similar to that described 
for the analogous conjecture in [psf 

As in Pq], there are two approaches to the asymptotics of E{n; s) for general n. We use the notation 



r{n; s) 



E{n; s) 



(1.25) 



In the first approach (see also |^ 27|) one differentiates (1.21) successively with respect to A. Using the known 
asymptotics of a{s;l) and the differential equation ( 1.21 ) satisfied by ct(s; A) for all A, one can find asymptotic 
expansions for the quantities 



O'n(s) 



5A^ 



and these in turn can be used to find expansions for the r(n;s). This approach is inherently incomplete since 
yet another undetermined constant enters the picture. And there are also computational problems since when one 
expresses the r(n, s) in terms of the (T„(s) a large amount of cancellation occurs, with the result that even the first-order 
asymptotics of r(n; s) are out of reach by this method when n is large. 
The second approach uses the easily-established identity 



r(n; s) = 



E 

'ii<...<i„ 



A,: 



(1-A.J---(1- A,J 



(1.26) 



where Aq > Ai > • • • are the eigenvalues of the integral operator K with A = 1 acting on (0, s). It turns out that this 
operator, rescaled so that it acts on (0, 1), commutes with the differential operator £ defined by 



Cfix) = (x(l 



^)/'(^))' - (^ + ?)/(^), 



with appropriate boundary conditions on /. Applying the WKB method to the equation, and a simple relationship 
between the eigenvalues of K (as functions of s) and its eigcnfunctions, we are able to derive the following asymptotic 
formula for the eigenvalues as s — *■ oo: 



1 - A,: 



27r 



T{a + i + l)i\ 



(1.27) 
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From this and (1.26) we deduce 



'■{n; s) - i Jl r(a + fc + 1) A:! I tt"" 2-"(2n+2a+i) ^-^-f » e^"^. 



(1.28) 



. fc=0 



For the special case a = 0, the quantity r(l;s) can be expressed exactly in terms of Bessel functions (see ( 2.30| ) 
below). 



II. DIFFERENTIAL EQUATIONS 
A. Derivation of the system of equations 



We shall use two representations for our kernel. The first is just our definition ( 1.2a ) using the notation ( |l.6| ) 



K{x,y) = 

The second is the integral representation 



x-y 



1 



This follows from the differentiation formula 
which gives the alternative representation 

^Ja + l{y/x)Ja{y/y) - Ja{y/x) y/y Ja+l{y/y) 



A- 



2 {x-y) 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



for K{x,y), and the Christoffcl-Darboux type formula (7.14.1(9)) of |^]. 
Our derivation will use, several times, the commutator identity 

[L,{I-K)-^] = (I-K)-^[L,K](I-K)-\ 

which holds for arbitrary operators K and L, and the differentiation formula 

^{I-K)-' - {I-Kr''^-^{I-K)-\ 
da da 

which holds for an arbitrary operator depending smoothly on a parameter a. We shall also use the notations 

M = multiplication by the independent variable, D = differentiation, 

and a subscript on an operator indicates the variable on which it acts. 

It will be convenient to think of our operator K as acting, not on J, but on (0, oo) and to have kernel 

K{x,y) Xjiy) 

where Xj is the characteristic function of J. Wc continue to denote the resolvent kern el of K by R{x, y) and note that 
it is smooth in x but discontinuous at y = aj. The quantities R{aj, aj) appearing in (1.5) are interpreted to mean 



hm R{aj,y), 



and similarly for pj and qj in formulas (1.7). The definitions (1.5) of u and v must be modified to read 



^Xj,{I-K)- 



)X^,(/-if)-V), 



(2.5) 
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where now the inner products are taken over (0,cx)). Notice that since 

(/ - Xy^i, = {I- K)-^^Xj in J 



for any function f, this agrees with the original definitions ( |l.8| ) of u and v. 
We have, by (p^), 

{{MD)^ + {MD)y)K{x,v) = t^^{<l>{xt)ct>{yt))dt ^-^<P{x)<P{y) - K{x,y). 

But it is easy to see that 

[MD, L] = {{MD)^ + {MD)y + /) L{x, y) (2.6) 
for any operator L with kernel L{x,y), where means "has kernel". Taking L{x,y) ~ K{x,y) Xj{y) gives 

[MD,K] = ^cP{x)^iy)Xj{y) - ^(-l)'=afci^(a:, a,) % - Ofe). 

(Recall the form ( |l.4| ) of J.) It follows from this and ( |2.3|) that 

[Mi?, (/ - = i Q(x) (/ - KY'Xj<I> (y) - Y^i-l^akRix, ak) piak,y), (2.7) 

where Q{x), and an analogous function P{x), are defined by 

Qix) := (/-i^)-V, ^'(a:) := (/-if)-V, (2.8) 

where p{x,y) — R{x,y) + 5{x — y) is the distributional kernel of (/ — i^)^^, and where A'* is the transpose of the 
operator K. (Note that K takes smooth functions to smooth functions while its transpose takes distributions to 
distributions.) Observe that 

qj = Q{aj), pj = P{aj). 



Next we consider commutators with M and use the first representation {2A) of K{x, y). We have immediately 

[M,K] = {c^[x)i:{y) - iAx)c^{y)) x/y). 



and so, by (2^) again, 

[M,(/-i^)-i] = Q{x){I-K')-^,pXj{y) - P{x){I-K'r^<(>Xj{y). 

Notice that since 

(/ - K'y^ijjXj = (/ - i^)" V = P on J, 

and similarly for Q, we deduce 



(2.9) 



In particular we have 



X — y 



R{aj,ak) = — (j ^ k) 

Gj — Qk 



R{x, x) = Q'{x) P{x) - P'{x) Q{x) (x g J). 



(2.10) 



(2.11) 



In order to compute R{aj,aj), and also the derivatives in (1.11) and ( 1.12| ), we must find Q'{x) and P'{x). We 
begin with the obvious 

X Q'{x) = MD{I - K)-^(j) {x) = {I - K)-^MD (j) {x) + [MD, {I - K)-^] (j) [x). 
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Using (2.7), and recalling (1.6) and ( p.5[ ), we find that 

xQ'{x) = P{x) + i Q{x) u - ^{-lfakR{x, a^) q^- 
Similarly, replacing by V' in this derivation gives 

xP\x) = {I- K)-HlD^j {x) + i Q{x) V - Y,i-^)'"akR{x, ak)pk. 
To evaluate the first term on the right side we use the fact that (p satisfies the difi'erential equation 



„2 J,// 



1 



(t)"{x) + X (j)'{x) + ^ (.T - a^) (t){x) = 0, 
which may be rewritten MDip (x) — ^(a^ — x) (j). Hence 

(/ - KyHlDi: (x) = ^ Q{x) - |(/ - K)-^M<t> {x) 



But we find, using (2^), that 

[M, (I - Xy^] 4> (x) = Q{x) V - P{x) u, 
and combining this with ( 2.13 ) and (|2.15| ) gives 

xP'{x) = i(a2 - x) Q{x) + i Q{x) v-^ P{x) u - Y^{-lfakR{x, ak)pk- 

It follows from ( p.U\ ), ( |2.12| ) and (|1|) that for a; G J 

aj R{aj,aj) ^ p^^ -j{a^ - +2v) g| + qj u + '^akR{aj,ak){qj Pk - Pj qk)- 



In view of ( ^.10 ) this is equation ( |l.l5| ) 



We now derive the differential equations (1.9)-(1.14). First, we have the easy fact that 

d 



dak 



K = {~lfK{x,ak)5{y-ak) 



and so by (2.4) 



dai 



-{I-K)-^ - {-lfR{x,ak)p{y,ak). 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



At this point we use the notations Q{x, a), P{ x, a) for P {x) and Q{x) to remind ourselves that they are functions of 
a as well as x. We deduce immediately from ( p. 17 ) and (2.8) that 



d d 

— Q(x, a) = {-l)''R{x, ak)qk, ^P{x, a) = {-l)''R{x, ak)pk- 
oak dak 



(2.18) 



Since qj — Q{aj, a) and pj — P{aj, a) this gives 



I* = (-l)'=i?(a„a,)g,, ^ = {-l)^R{a,,ak)pk, {j ^ k). 
oak oak 



In view of ( pl0| ) again, these are equations (|1.9|) and ([1.10[). Moreover 



dqj _ ^ d ^ d -^Q(^^ 
daj dx daj ' 



dp 



d d 



da 



'dx da 



Equations ( pi) ) and ( p^ follow from this, ( |2l8| ), ( pl^ ) and ( plc| ). 

Finally, using the definition of u in (2^), the fact 



and (2.17) we find that 



But 



{(f>Xj,R{-,aj))= J R{x, aj)(l){x) dx = J R{aj,x)(l){x) dx 
since R{x, y) — R{y, x) for x, y G J. Since R{y, a;) = for x ^ J the last integral equals 

R{aj,x)(j){x) dx = Qj — 4>{cLj)- 



This gives (1.13), and (1.14) is completely analogous. 

We end this section with two relationships (analogues of (2.18) and (2.19) of |^^) which would allow us to express 
u and V in terms of the qj and pj if we wished to do so. (They will also be needed in the next section.) These are 



(2.19) 



= II(-l)' - (a' - aj + 2v)q] + 2pjq,u) 



(2.20) 



To obtain the first of these observe that (|l.9| ) and ( |l.ll| ) imply 



E 



d 
dak 



1j = P] + in^^ 



while from (1.13) and (1.14), 



^ (2i; + -y) = 2(-l)^gj(p, + -a^u). 



If we multiply both sides of the previous formula by (— l)^ajgj and sum over j what we obtain may be written 



or equivalently 



_d_ 

da 



5 \ 1 2 

OfcT; — (2w + -M + 7i . 



It follows that the two sides of ( ^.19| ) differ by a function of (ai, . . . , a2m) which is invariant under scalar multiplication. 

Since, as is ea sily seen, both sid es van ish when all aj — their difference must vanish identically. 

To deduce ( ^.2C| ) we multiply ( 2.17 ) by ak and sum over k and then add the result to ( ^.7| ), recalling ( |2.6| ), to obtain 
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for x,y £ J. This gives 



9u 
9a,, 



If we multiply both sides of this by (—1)-' and sum over j we deduce, by an argument similar to one just used, that 

^(-l)%,i?(a„a,) = iw. (2.21) 



Substituting for ajR{aj,aj) here the right side of ( 1.15 ) we see that the resulting double sum vanishes, and ( 2.20| ) 
results. 



B. The ordinary differential equation 



In this section we spe cializ e to the case J = (0, s) and derive (among other things) the differential equation ( 1.16| ) 
and the representation ( 1.19| ). In the notat ion o f the la st section m = 1, ai = 0, 02 = s. We shall write q{s),p{s), R{s) 
for q2,p2, R{s, s), respectively. Equations (1.11) (1.14) become 



sq = p + -qu 

, 1 . 2 ^ 1 1 
sp = -[a -s)q+-qv--pu 

u' = q^ 
v' — pq. 



It is immediate from (2.21) and (2.24) that 



{sR{s)y = -q{sf, 



(2.22) 

(2.23) 

(2.24) 
(2.25) 



(2.26) 



and since 



ds 



log D{J;X) = -R{s) 



(see (1^)), we obtain the representation ( 1.1£ ) 



To obtain the differential equation (1.16) we apply to both sides of (2.22) and use (2.22), (2.22) and (2.24) 
What results is 



s {s q'Y ^\{a^~ s)q + + 8v)q + q\ 



(2.27) 



But ( 2.19 ) in this case is 



and so the above can be written 



+ 8u = 4s - 4u, 



s{sq'y = jia^ - s)q~ ^uq+^sq^. 



(2.28) 



Next, we square both sides of (2.22) and use ( 2.20] ), which now says 

u = 4p^ - (a^ - s + 2-y)g^ + 2pqu, 

and find that 

isqr^\u + \{a'-s)q' + ^^q'{u' + 8v). 
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Combining this with ( 2.2l\ ) gives 



and combining this with (2.28) gives the desired equation (1.16). The boundary condition (1.17) fohows from the 
Neumann exp ansion of the defining expression (1.7) for q. 
Using (1.18) one easily verifies that 



1 



■ cos ^|: + s 



ds 



, csc^ i/" + - 
4s s 



-cos^V + s -r- cot^V+ C-— - 

4 \ds 4s 4 s 



satisfies ( 2.26| ) where c is a constant of integration. That this constant is equal to a^/4 follow s from the small s 
expansion of R(s). (U se th e fact that for s — > 0, R{ s) ^ K(s, s) and that, as follows from ( 1.2b| ), there is no simple 
pole in s.) Equation ( 1.21 ) follow s from ( 1.20 ) and ( 1.18|) . 

Here is the argument why (1.16) must be reducible to some Painleve equation (or one of the other simpler differe ntia l 
equa tion s on Painleve's list). The derivation of ( 1.16 ) used only the facts that the Bessel ker nel h ad both forms (2.1) 
and (2.2) and that the function cf) satisfied the differential equation (2.14). (Of course tp in ( ^.l[ ) must be defined as 
MD(f).) This equation has a 2-complex-parameter family solutions and this gives a 2-coniplex-parameter family of 
kernels defined by (2.1). They can be shown to satisfy (2.2). We replace the kernels K{x,y) by sK{sx,sy) and have 
them act on (0, 1) rather than (0, s). These operators on (0, 1) depend analytically on the complex variable s (except 
for a branch point at s = 0) and the corresponding g(s) can have, aside from a branch point at s = 0, only poles 
which occur at the values of s for which A = 1 is an eigenvalue of the operator. (The resolvent of an analytic family of 
compact operators has a pol e whenever A = 1 is an eigenvalue.) Thus the general solution (i.e., 2-complex-parameter 
family of solutions) of ( 1.16| ) has only poles as moveable singularities. Since the equation is of the form q" = rational 
function of q' and g, it must be reducible to one of the Painleve types. 

We mention that this argument requires |q;| < 1 since it is only then that all solutions of (2.14) give compact, or 
even bounded, operators on L2- For other a it may be that we just have to replace L2 by an appropriate space of 
distributions. 



C. r(l;s) for a = 



If we set a = and make the change of variables s ~ x"^ , the differential equation for %j: (recall (1.18)) becomes 



^" + i^' = 1 sin(2V') 
X 2 



(2.29) 



and we want the solutions holomorphic at the origin. The linearization of this differential equation is the modified 
Bessel equation and all solutions of the linear equation are linear combinations of Iq{x) and Kq{x). Flaschka and 
Newell |1(| have shown, using metho ds of monodromy preserving deformations and singular integral equations, that 
the general 2-parameter solution to ( 2.29 ) can be viewed as a "perturbation" of this linear solution. (Precisely, they 
derive a singular integral equation whose Neumann expansion in a particular limit gives ip{x) — the fir st te rm in this 
expansion is a linear combination of Bessel functions.) The one-parameter family of solutions to ( 2.29| ) that are 
holomorphic at the origin has the representation |10 



V'(x;/i) = fiTpiix) + -^^'3(2^) + -^^5(2;) + • • • 

where ijji{x) = Iq (x ) , and /i^ = 1 — A. Note that we are using the slightly confusing notation ip^x; ^) to denote the 
function t/j of (I.IS) after the change of variables s = x'^ and /i^ = 1 — A. Multiple integral representations (obtained 
from a Neumann expansion) for the higher 'ipj's can be easily derived from | [To[ |. 
The resolvent kernel i?(s), s = x'^, is given in this special case by 



4i?(s) 



COS' 



= 1 + (1 - A) {i4>[y - i^!) + (1 - A)^ (^4 - Ms + ^'iV'a) + 
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Thus (recall (0|) and (lOsI)) 



where 



Therefore 



r(l;s) 



Ri{t)dt 



d 



R^it) = —R{t) 



r(l;s) 



1 
1 



x{ijl{x)-{^[{x)f) dx 



X {lo{xf - ll{x)) dx 



(2.30) 



The last equahty follows from (5.542) of |12|. 

We point out the curious fact that (after letting ijj — + iij)) the same differential equation ( ^.29| ) and closely related 
T- function arise in the 2D Ising model p0|,pi|,p6[ except that here the boundary condition is il}{x) ~ iJ.Ko{x) as a; — > c». 



III. ASYMPTOTICS 



A. Asymptotics of the cr-equation 

In the case of the finite iV ensemble and a = 0, Edelman and Forrester ||ll[] (by a direct evaluation of the 
integrals defining the probability En{0; s)) has shown that ^^(0; s) is exactly equal to e~^/"'. From this it follows, for 
a — and A = 1, that a{s; 1) — s/4. For general a and A = 1 it is therefore reasonable to assume an asymptotic 
expansion of the form: 



a{s; 1) = Cis + C2S^^^ + cs + c^s 



-1/2 



oo. 



Substituting this into the differential equation (1.21) results in equations that uniquely determine the coefficients Cj 
once a choice in the sq uare root Vc^ is made. Since for a = our asymptotic expansion of E{0; s) must agree with 
those of Dyson (recall (|l.3|)), we see that we must choose the square root —a —a. A calculation then gives 



a{s;l) 



4-r 



1/2, 



a a 1/9 a 
-+16^ 



2048 



(128a'' + 720a^ 



2 

' 16 
225)s"^/2 



1 , ^ /ic„2 , A\„-3/2 I ^/^„.2 , n^„-2 



^(16a^ + 9).- 



+ _(4a^ + 9)s- 



(3.1) 



from which (1.24) follows. 



B. Asymptotics via the commuting differential operator 

Throughout this section we take A = 1. The operator K, when rescaled to act on (0, 1) instead of (0, s), has kernel 
s K{sx, sy). By (^) this is equal to 

1 

4>{sxt) 4>{syt) dt, 
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and so K (rescaled, as it will be throughout this section), is the square of the operator on (0, 1) with kernel 

Jix,y) = —(f>{sxy) = —Jaiy/sxy). 



This function satisfies the differential equation 



2 / sxy cx, , 

X J^oo + xJx + —)J = 0. 

The operator J will commute with a differential operator 

C ^ ^a{x)^ + P{x) 
ax ax 

if a(0) ail) = and if 

a{y)Jyv + a'{y)Jy + j3{y)J = a{x)Jxx + a'{x)J^ + (3{x)J. 
If we use (^.2[) we see that this will be satisfied if 



-a{x) 



+ a'{x)Jx + /3(x)J 



(3.2) 



= the same expression with x and y interchanged. 
Equating the terms involving the first derivatives of J gives 

(a' (a;) - x^'^a{x))Jx = [a' {y) - y~'^a[y))Jy. 
But X Jx = y Jy, so the above will hold if 

x'^(a'{x) - x^^a{x)) = y^'^{a'{y) - y~'^a{y)). 
This is satisfied if a[x) is a quadratic without constant term and, of course, we choose 

a{x) = x{l — x). 

What is required of /3, then, is seen to be 



which is satisfied by 



m - x{i - .)(^ ~^)^ m - yii - y)(g ~ 



a" sx 



We write the differential equation, for which the eigenfunctions f{x) are the eigenfunctions of J, as 

(x(l - x)fx)y + i^^s-^- ■^)f{x) = 0. (3.3) 

The boundary conditions are that f{x) be bounded as a; ^ 1 and that f{x) be asymptotic to a constant times x"^"^ 
as X —^ 0. The reason we wrote the eigenvalues as we did is that for each i the ^ corresponding to the i'th largest 
eigenvalue is bounded as s ^ oo. This is easily seen by an oscillation argument. So we assume i is fixed and proceed 
to find the asymptotics of the corresponding cigcnfunction f{x) as s ^ oo. We assume it normalized so that 

fix) - x"/^ as a; -> 0. (3.4) 
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1. The region a; ^ 1 



The approximating equation is 



The solution of this equation which satisfies (|3.4[) is 



X Ax 4 



where $ is the confluent hypergcomctric function H. We deduce that when a; <C 1, 



f{x) ~ a;"/2 V^-/2 - /x, 1 + a, V^x). (3.5) 



2. The region a; ^ 1, ^/sx ^ 1 

Case 1. ^ ^ i + (i = 0, 1, . . .). Then from the known asymptotics of $ as its argument tends to oo (6.13(3) of 
[^) we deduce that 

/(a;)^ J^(i±^s-^(^+'^)a;-^-^eV^-/2. (3.6) 



Case 2. /i Then 



<f>(-z, 1 + a, ^sx) = ^> J^. ^ Lti^sx) 
1 (a + I + 1) 



(3.7) 



where Lf is the generalized Laguerre polynomial (6.9(36) of ^). So we find that in this case 

fix) - (-l)^-F?^^^ x'^^ e"^"/'- (3-8) 
r(Q; + « + 1) 



3. The region a/s a; » 1, (1 — a;) s > 1 

Here we use the standard WKB approximation for the solutions of a differential equation {p f')' + qf = given 
by /{pq)i where y = J \fqjv- In the case of our equation ( |3.3| ), the range of validity of the approximation is as 
indicated in the heading of this section. To be definite, we take 



y = - y/q{z)/p{z)dz. 

J X 

It is easy to compute that for the range in question we have 



y = ~Vs y/l — X ~ fi log a; + 2 ^log(l + \/l — x) + o(l). 
Also, in this case is asymptotically — s/4 times x^{l — x). Hence 

fix) - (a(s) x-f" (1 + 71^)2*^ ^-V^VT^ 

+b{s) x" (1 + VT^y^" e^'^^ (3.9) 
where a(s) and b{s) are constants depending on s. 
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4- The region x 1 

Letting y — 1 — f{x) = g{y) gives the approximating equation 

{ygiv)')' ~'-^9{y) = o. 

The general solution of this is a constant times Io{-\/y) where /q is the modified Bessel function. Thus we deduce that 
as X — > 1 

f{x) ^ c{s) lo (Vs(l-a^)) (3.10) 

for some c(s). 



5. Determination ofa{s), b{s) and c{s) 



From the asymptotics of Iq at infinity and ( 3.10 ) we deduce that 

c(s) ev^^ 



2tt Si (1 — a;) 4 



when a; — > 1 and s(l — x) — > oo. Comparing this with (3.9) shows that a{s) ~ and that 



c(s) ^ V2t: s-i b{s). 



(3.11) 



And now comparing (3^) with (3^) and (3^) in their overlapping range of validity we see that we must be in case 2 
and that 



bis) ^ i-iy J^"" + , 22"+"+i e-^. 
^ ^ ^ ' T{a + i + l) 



(3.12) 



Wh at we shall need from all this is, first, the asymptotics of /(I). This follows immediately from ( 3.10| ), ( |3.11 ), 
( |3.12D , and the fact that I„{0) = 1: 



/(I) ^ i-irV2^-£iq±^ si+ie-^2'"+'^+\ 
r{a + 1 + 1) 



(3.13) 



We shall also need the asymptotics of / f{x)'^dx. It follows from the asymptotics we have derived that the main 
contribution to this integral comes from an arbitrarily small neighborhood of x = 0. It follows from (3.5), (3.7) and 
the fact 



x"e^^ Lf{x)^dx 



T{a + i + l) 



that 



fixfdx 



r(a+l)2i! 



S 2 2 . 



r(a + z + 1 

We put these two relations together to get what we really want, which is 



/g/(x)2dx r{a + i + l)i\ 



i+f +1 „-2^/s r,4i+2a+2 



(3.14) 
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6. The asymptotics of \i 



Since C commutes with our integral operator K, rescaled to act on (0, 1), the set {/;} of eigenfunctions corresponding 
to the eigenvalues /io < Mi < ... of £ is the set of eigenfunctions of K corresponding to its eigenvalues, in some order. 

Lemma 1. The eigenvalues of K are simple. 

Proof. Since the eigenvalues of K are the squares of the eigenvalues of J what we have to show is that if /i (x) and 
f2{x) are eigenfunctions of C (they need not correspond to any particular /i here) and if for some v we have either 

1 

Ja{y/sxy)Ji{y)dy ^ vfi{x) (i = l,2) 





Ja{\/sxy)fi{y)dy ^ vfi{x), / Ja{\/sxy)f2{y)dy ^ -v f2{x) 

"'0 

then /i and /2 are linearly dependent. We shall assume them normalized so that they both have the value 1 at a: = 1. 
Notice first that v ^ Q for otherwise if we expand the /i equation near a; = we would find that 

1 

+k 



fi{x)dx ^ Q, fc = 0, 1,. 







and /i would be identically 0. 
Next we define 



H{x) = e-^/2j„(V^), g^{x) = e-^/V.(e""), 
make the obvious variable changes, and find that our relations become 

H{x + y)g,{y)dy^vg,{x), (z = l,2), (3.15) 



H{x + y)gi{y)dy ^ V gi{x), / H{x + y) g2{y) dy ^ ~vg2{x). 



(3.16) 



What comes now is almost identical to the proof of Lemma 1 of |g8|. Assuming first that (3.15) holds, we differentiate 
twice this relation with i = 1 and then integrate by parts twice to obtain 

ygi{x) = -H'{x)+H{x)g[[0)+ ( H{x + y) g,{y) dy. 



If we multiply both sides by 32(2;) and integrate, using ( 3.15| ) and its differentiated version, we obtain (recall that 

9^(.0) = 1) 

V g';{x)g2{x)dx^-vg'M + vg[{Q) + v g2{y) g'Hy) dy. 



Thus, since v 0, we have 5^(0) = g2{0)- Equivalently, /{(I) = /2(1)- But since also /i(l) = /2(1), it follows from 
equation (2) that the corresponding eigenvalues /i must be the same, and so the eigenfunctions are the same. 
Next, assume ( 3.16| ) holds. Differentiating both sides of the first relation once and integrating by parts give 



yg[{x) = -H{x)- H {x + y) g[{y) dy . 
Jo 

Multiplying both sides of this by 32(2;) and integrating, using the second part of ( 3.16|) , we obtain 

/■oo poo 

^/ g'i{x) g2{x) dx = -V + V I g2{y)g'i{y)dy, 
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contradicting i/ 7^ 0. 



Now that we know the eigenvalues Xi of K are simple we can order them so that Aq > Ai > . . . . There is a 
permutation cr of N such that the eigenvalue corresponding to fi is Xa{i) ■ This permutation is independent of s since 
the eigenvalues fii and A^ both vary continuously with s. The next lemma allows us to compute dXi/ds for each i. 

Lemma 2. Let X be a simple eigenvalue of an integral operator on (0, 1) with symmetric kernel of the form s K{sx, sy) 
where K{x, y) is smooth and symmetric. Let f he a corresponding eigenfunction. Then 

dX X f{lf 



Proof. We may assume / normalized so that f{x)^dx — 1. Then according to Lemma 3 of we have 

~r = f i ^{sK{sx,sy)) f{y) f{x)dydx. 
ds Jo Jo OS 

(This holds no matter what the form of the kernel.) In the case at hand the first factor in the integrand equals 

d d 
K{sx,sy) + {x— + y—)K{sx,sy) 

and we easily deduce from this that the integral itself equals 



( [ f{xfdx+ [ xf{x)f{x)dx+ [ yf'iy)f{y)dy 
\Ja Jo Jo 



and integration by parts shows that the expression in parentheses equals /(l)^. H 
Remark. It is easily seen that the conclusion of the lemma holds when the kernel has a mild singularity at such as 
our Besssel kernel has when a < 0. 

Lemma 3. For each i we have Xi ^ 1 as s ^ 00. 

Proof. The Hankel transform, when rescaled by the variable change x 1-^ y/x, is the integral operator H on (0, 00) with 
kernel ^Ja{y/xy) and so our operator K on (0, s) may be thought of as Pg H PsH Ps where Ps denotes the projection 
from L2(0,oo) to ^2(0, s). Since, as is well-known, = /, the minimax characterization of the eigenvalues shows 
that for each i the eigenvalue of Ps H Ps with i'th largest absolute value tends to ± 1 as s — > cx). Since our operator 
is the square of this one, the statement of the lemma follows. ■ 

We can now deduce the asymptotic formula (15) for th e eig envalues A^. We apply Lemma 2 to the eigenvalue Acr(i) 



associated with the eigenfunction fi{x) of (p.3|) and use (3.14) to deduce 



ds T{a + i + l)i\ 

Recalling Lemma 3 we see that we can integrate from s to 00 and we obtain 



"^'^ r{a + i + l)i\ 

It remains to show that a{i) = i for all i. But it is clear from the above formula that i < j implies that Xa(i) > A(j(j) 
for large s (and so for all s) and therefore a{i) < cr{j). Since tj : N ^ N is onto, we must have a{i) = i for all i. 
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